Slow Excitation Trapping in Quantum Transport with Long-Range Interactions 
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Long-range interactions slow down the excitation trapping in quantum transport processes on a one- 
dimensional chain with traps at both ends. This is counter intuitive and in contrast to the corresponding classical 
processes with long-range interactions, which lead to faster excitation trapping. We give a pertubation theoreti- 
cal explanation of this effect. 
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Building a quantum system from scratch has become pos- 
sible due to recent experimental advances in controlling and 
manipulating atoms and molecules. It has actually become 
possible to tailor theoreticians favourite one-dimensional sys- 
tems using, e.g., ultra-cold atoms in optical lattices, see ^ 
and references therein. From a dynamical point of view, this 
allows for these systems to compare the theoretical predici- 
tons for the transport of charge, mass, or energy to the ex- 
perimental results. In turn, the experimental findings might 
eventually lead to a refinement of the theoretical models. 

The tight-binding approximation for the transport of a 
quantum particle over a regular structure (network) is a sim- 
ple description which is equivalent to the so-called continous- 
time quantum walks (CTQW) with nearest-neighbor interac- 
tions (NNI) 12, [It] . Recently, several experiments have been 
proposed addressing CTQW, e.g., based on wave guide arrays 
p], atoms in optical lattices jsl |6^, or structured clouds of 
ultra-cold Rydberg atoms In some of these experiments 
one finds long-range interactions (LRI), such as in Rydberg 
gases, where also blockade |9] and antiblockade |10] effects 
have to be considered. In a recent study of the effect of LRI 
on the quantum dynamics in a linear system it has been found 
that CTQW for all interactions decaying as i?"" (where R 
is the distance between two nodes of the network) belong 
to the same universality class for v > 2, while for classi- 
cal continuous-time random walks (CTRW) universality only 
holds for > 3 H. 

Coupling a system to an absorbing site, i.e., to a trap, al- 
lows to monitor the transport by observing the decay of the 
survival probability of the moving entity, say, the excitation. 
In the long-time limit and for NNI the decay is practically ex- 
ponential for both, classical systems modeled by CTRW ifllll 
and quantum systems modeled by CTQW ||3.[l2l- At interme- 
diate times, which are experimentally relevant, there appear 
considerable, characteristic differences between the classical 
and the quantum situations |7]. 

Here, we study the quantum dynamics of one-dimensional 
CTQW with LRI in the presence of traps and use the similarity 
to CTRW for a comparison to the respective classical case. 
Without traps, we model the quantum dynamics on a network 
of connected nodes by a tight binding Hamiltonian Hq. For 
the coiTesponding classical process, we identify the CTRW 
transfer matrix Tq with Ho, i.e.. Ho = — Tq; see e.g. [2, 
[stl for details. For undirected networks. Ho is related to the 



connectivity matrix Ao of the network by Ho = Ao. When 
the interactions between two nodes go as R^'^, with i? = |fc — 
j\ > 1 being the distance between two nodes j and k, the 
Hamiltonian has the following structure: 
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We restrict ourselves to extensive cases > 1), i.e., we ex- 
plicitly exclude ultra-long range interactions. The correspond- 
ing NNI Hamiltonian is obtained for = cx), in which case 
only the leading terms with i? = 1 do not vanish. 

The states \j) associated with excitations localized at the 
nodes j {j = 1, • ■ • , N) form a complete, orthonormal basis 
set of the whole accessible Hilbert space ({k\j) — dkj and 
Efe 1^ )( ^1 ~ l)- In general, the transition probabilities from 
a state \ j) at time to = to a state |fc) at time t read Hkjit) = 
\akj{t)\^ = exp[— iHo(i^)t]|j)|^. In the corresponding 
classical CTRW case the transition probabilities follow from 
a master equation as pfcj (t) = (fc| exp(Tot)|j) iQIsll- 

Now, let the nodes m (m e M and M C {!,..., N}) be 
traps for the excitation. Within a phenomenological approach, 
the new Hamiltonian is H(z^) = Ho(i^) — iF, with the trapping 
operator iT = iV J2meM 1"^) 01 for details. As 

a result, H is non-hermitian and has N complex eigenvalues. 
El = d - i^i (I = 1, . . . , N) with 7; > 0, and N left and N 
right eigenstates, denoted by j^I';} and respectively. The 
transition probabilities follow as 



TTfcjW = |Xlexp(-7ii)exp(-ieit)(fc|*i)(^'i|j) 



(2) 



where the imaginary parts ji of Ei determine the temporal 
decay. For the incoherent classical process the description 
by CTRW is quite similar: The new transfer operator reads 
T(i^) = Toil/) - r = -Ao(i^) - r, which is real and sym- 
metric, leading to the eigenvalues —A; (A; > 0) and corre- 
sponding eigenstates |$;). Note that due to the different in- 
corporation of the trapping operator in T(i') and H(i^) the 
coiTesponding eigenvalues and eigenstates will differ With- 
out trapping we have To(j^) = — Ho(i') and thus A; = Ei and 
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In order to make a global statement for the whole network, 
we calculate the mean survival probability for a total number 
of M trap nodes, 



1 



N - M 



(3) 



i.e., the average of TTkjit) over all initial nodes j and all final 
nodes k, neither of them being a trap node. Classically, we 

will consider Pa/ W = - M) E,^a^ Efe^A^ Pfej W- 

For intermediate and long times and a small number of trap 
nodes, IlM{t) is mainly a sum of exponentially decaying 
terms 101: 
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1=1 



exp(-27,i). 



(4) 



If the imaginary parts 7; obey a power-law with an exponent /i 
(7; ^ al'^), the mean survival probability scales as nj\/(t) ^ 
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FIG. 1; (Color online) !/-dependence of the quantum mechanical 
Umit) and the classical PM(t) decay behaviors for a chain of 
N = 100 sites; here (a) T = 0.001 and (b) T = 1. The inset in 
(a) shows a close-up picture of the region where IlMit) and PMit) 
cross. The inset in (b) shows power-law fits to IlM(t) in the interme- 
diate time regime with exponents l//i, where the /i are taken from 
Fig.Hb). 

In Ref. ifT^l an experimental setup was proposed, which is 
based on a finite linear chain of clouds of ultracold Ryd- 
berg atoms with trapping states at both ends (m = l,iV). 
There, the dynamics was approximated by a NNI tight bind- 
ing model, which - for a ring without traps - has been shown 
to behave in the same fashion as systems with LRI of the form 
R^'^ for which > 2 [8J. The Rydberg atoms interact via 



dipole-dipole forces, i.e., the potential between two atoms de- 
cays roughly as R^^. 

For the finite chain with m — 1,N, Fig.[T]shows a compar- 
ison of the quantum mechanical Il]\j{t) and of the classical 
Pm (t) behaviors for different i' and F, which were obtained 
by numerically diagonalizing the corresponding Hamiltonian 
H(i^) and transfer matrix T(zy), respectively. Clearly, for both 
F-values the LRI lead to a slower decay of IlM{t), i.e., to a 
slower trapping of the excitation, which is counter intuitive 
since the opposite effect is observable for classical systems 
where the decay of PM{t) becomes faster for decreasing ly, 
see below. By increasing the trapping strength F, the differ- 
ence between the quantum and the classical behavior becomes 
even more pronounced, compare Figs. [Tfa) andflfb). Gener- 
ally, for IIm (t) the change in F results mainly in a rescaled 
time axis, since the imaginary parts 7; are of the same order 
of magnitude when rescaled by F. For the specific case of the 
Rydberg atoms (i^ = 3 and F = 1) one observes the largest 
difference between the Unlit) and the PAi{t) behaviors. To 
understand this phenomenon, we continue to analyze Hm (t) 
within a perturbation theoretical treatment. 

When the strength of the trap, F, is small compared to the 
couplings between neighboring nodes, we can evaluate the 
eigenvalues using perturbation theory, see, for instance. 111 311 . 



Let 1*1°''} be the hh eigenstate and E^^^ G M be the /th eigen- 
value of the unperturbed system with Hamiltonian Ho {i^). Up 
to first-order the eigenvalues of the perturbed system are given 
by 



(0) 



(5) 



Therefore, the correction term determines the imaginary parts 
7;, while the unperturbed eigenvalues are the real parts e; = 
-E; . Having only a few trap nodes, the sum in Eq. ^ con- 
tains only few terms. Moreover, from Eq. (|5]l we also see that 
the imaginary parts 7; are essentially determined by the eigen- 
states 1*1"^} of the system without traps. A change in these 
states will also lead to a change in the 7;. As we proceed to 
show, this is exactly what happens by going from NNI to LRI. 

Without loss of generality, an eigenstate of a finite chain 
with NNI can be written as (/ = 1, . . . , TV) 



1 ^ 



I = N 



N 



(6) 



^cos[(2j-l)0,/2]|j) else, 



where for convenience we take 9i = tt{N — G [0, 7r[; 

the corresponding eigenvalues are Ef'^ =2 — 2 cos 9i (note 
that the smallest eigenvalue is E^^^ = 0). Thus, to first order 
perturbation theory we obtain from Eqs. (|5]l and (|6]l as imag- 
inary parts 7JV = 2F/iV and 7; = IaT/N)cos^ (61/2) = 
(2F/A^) [1 + cos 61;] for ? = 1, TV - 1, which for / < iV 
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yields 7/ ~ Z^. In this case the mean survival probability will 
scale in the corresponding time interval as 11^(1) t^^^^. 

Formally, we can perform the continuum limit N ^ 00 (by 
taking now AT/N = a finite). Then the sum in Eq. ^ turns 
into an integral such that 



(a) 



nA/(<) 







e-'^'loiat), (7) 



where Iq (at) is the modified Bessel function of the first kind 
1 1411 . From this we get for large t that nA/(t) ~ t^^^^, which 
confirms the previous results. Note, however, that for small N 
the smallest 7; -value is finite and, therefore, the scaling of 7/ 
holds only in a quite small interval of ^-values. Hence, also the 
time interval in which Umit) scales with the exponent —1/2 
is rather small. A lower bound for scaling is given by the 
behavior of 7/ for I w N/2 (corresponding to smaller times 
than for I <C N). Here, 7/ is linear in I, which leads to a 
lower bound of fi > 1 for the scaling exponent. An exponent 
fi which is valid over a larger ^-interval will therefore be in the 
interval [1, 2] and, consequently, the exponent for IlM{t) will 
lie in the interval [—1,-1/2]. 

In the case of periodic boundary conditions, one finds 
translation-invariant Bloch eigenstates regardless of the range 
of the interaction fs^l. In the case of Eq. ([U, however, the 
eigenstates for LRI differ from the ones for NNI [Eq. ©J; In 
Eq. ([T]i the finite extension of the chain destroys the transla- 
tional invariance. As is immediately clear from Eq. (|5]), this 
also implies that the imaginary parts of the eigenvalues, eval- 
uated based on first order perturbation theory, will change. 

For large exponents i' we can regard the LRI as a small per- 
turbation to the NNI, i.e., having Ho(j^) = Hq + H^, where 
Hi, contains only the correction terms to the NNI case Hq. 
This allows us to calculate from the unperturbed states | ) 
the perturbed eigenstates j'!';) up to first order. Taking the 
states I 'J;) to be the eigenstates of the LRI system without 
traps, we readily obtain the imaginary parts 7; for small trap- 



ping strength from Eq. (0) as 7; = 2r| 



,(0)\ 



', where 



It is straightforward, although cumbersome, to calculate 
the corrections to the imaginary parts 7; from Eq. (O. For 
large v the coupling to the next-next-nearest neighbor is by 
a factor of (3/2)'' smaller, for v — 10 this is about one 
and a half orders of magnitude. Taking, for fixed v, only 
nearest and next-nearest neighbor couplings into account al- 
lows us to obtain simple analytic expressions. The pertur- 
bation term H,, is now tri-diagonal. Its non-zero elements 
are (j - 2|H,|j) = {j + 2|H,|j) = -2"'' and its diago- 
nal elements follow from (j|H,y|j) = — thus 
(jIH.Ij) = 2- for 2 < J < iV - 1 and {3\RA]) = 2-+i 
else. We hence obtain from Eq. (|8]l 
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FIG. 2: (Color online) Correction term — 
100 and for (a) = 10 and (b) 1/ = ^. The direct numerical eval- 
uation (solid black line) is compared to the perturbation theory ex- 
pression Eq. ([8} (dashed-dotted green line) and to the approximate 
expression Eq. l|9ll (dashed red line). 

Figure |2] shows the difference (1|*/) - (1|*|°^) for N = 
100 and for (a) v = \Q and (b) v — ^. The numerical exact 
value (solid black line) is obtained by computing separately 
(1|5';) and and subsequently taking the difference; 

the result is then confronted to Eq. (O (dashed-dotted green 
line), determined numerically, and to Eq. (|9]l (dashed red line). 
For V = 10, the agreement between all three curves is remark- 
ably good, see Fig. I2a), which justifies the assumptions lead- 
ing to Eq. (|9]l. For smaller [z^ = 5 in Fig.|2b)] there is still 
a reasonable agreement between Eq. (O and the exact result; 
however, taking only nearest and next-nearest neighbors into 
account leads to evident deviations, see the dashed red line in 
Fig.Hb). 

Now, from Eq. (|9]l we get 



^ (0) , 
11 « 7i + 



+ 0(2-2"), 



(10) 



where 7^°^ is the NNI expression given above and 7|^-' — 
[ST/N) cos {61/2) sin (26*;) sin (6';/2) the coiTection due to 
the LRI. Again, the smallest 7; -values are those for which I ^ 
N, which leads to a decrease of the imaginary parts 7; because 
7p^ < for I ^ N . Here, one can approximate the imagi- 
nary parts by a power-law, i.e., 7/ ~ Z^. A rough estimate of 
the scaling exponent /i, assuming ^ 1 can be readily given. 
For this we note that from Eq. (fTOl i we have In 7;+i — In 7; « 



More- 



over, the term = [\iii[% - ln7;^°^] / [ ln(/ + 1) - In/] 
gives the exponent for the NNI case and the term /i'^^ = 

b\+ih\+i - l\'^h\'^]l[Hl + 1) - In/] is the LRI cor- 
rection. Thus 



ln7;-|-i - In 7; 
ln(/ + l)-ln/ 



(0) 



(1) 



(11) 



Since /i^^^ is strictly positive for small /, the inclusion of LRI 
leads to a decrease of 7; when compared to the NNI case. In 
turn, this results in a slower decay of I{M{t)- 
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FIG. 3: (Color online) Imaginary parts 7; (dots) in ascending order 
for LRI systems with u = 2,S, 4, and for NNI for = 100 and (a) 
r = 0.001 and (b) T = 1. 



Figure [3] shows the imaginary parts 7; for a chain of = 
100 nodes with LRI (v = 3, 4, 5) and with NNI. For small I 
and NNI, the 7; obey scaling with the exponent /i = 2, as dis- 
cussed above. Introducting LRI, i.e., decreasing v, increases 
the scaling exponent to /i > 2. Consequently, the scaling ex- 
ponent l//i for n7\/ (t) decreases, leading to a slowing-down 
of the excitation trapping due to LRI. 

In the classical case decreasing v leads to a faster excita- 
tion trapping, which is obervable in a quicker decay of Pnj{t). 
This can also be deduced from a perturbation theoretical treat- 
ment. As can be seen from Fig.[T](see also Fig. 2 of Ref. 101), 
the decay of PM{t) is exponential already at intermediate 
times and is dominated by the smallest eigenvalue Ajv and the 
corresponding eigenstate |$Ar) of the transfer operator T(zy): 



N 



PM{t) - ^^r^Eexp(-A,i)| Y.m'^i) 

1=1 k^M 



(12) 



Calculating A^v and the prefactor | J2k^M (M^n) \ for large 
ly and small T shows that with decreasing v the smallest 
eigenvalue Xn increases while the prefactor decreases. To- 
gether, this confirms our numerical result of a quicker decay 
for PmW, see Fig. □ 

Finally, we comment on the impact of our results on the ex- 
periment proposed in Ref. |i7J . Here, clouds of laser-cooled 
ground state atoms are assembled in a chain by optical dipole 
traps which are then excited into a Rydberg S-state, see 



10] for details. The Rydberg atoms interact via long-range 
dipole-dipole forces which is advantageous in many ways. As 
can be deduced from Fig.[T] the time intervals over which the 
decay follows the power-law are enlarged by the LRI. For 
ly = 3 the transition to the long-time exponential decay oc- 
curs at times which are about two order of magnitude larger 
than the ones found for the NNI case. The difference between 
a purely coherent (CTQW) and a purely incoherent (CTRW) 
process is enlarged due to the LRI, allowing for a better dis- 
crimination between the two when clarifying the nature of the 
energy transfer dynamics in ultra-cold Rydberg gases. 

In conclusion, we have considered the quantum dynamics 
of excitations with LRI on a network in the presence of ab- 
sorbing sites (traps). The LRI lead to a slowing-down of the 
decay of the average survival probability, which is counter in- 
tuitive since for the corresponding classical process one ob- 
serves a speed-up of the decay. Using pertubation theory ar- 
guments we were able to identify the reason for this slowing- 
down; it results from changes in the imaginary parts of the 
spectrum of the Hamiltonian. 
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